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CN ' Abstract 



The possibility of dynamical stabilization of an internal space is investigated for a multidimensional 
cosmological model with minimal coupled scalar field as inflaton. It is shown that a successful dynamical 
compactification crucially depends on the type of interaction between the geometrical modulus field and 
the inflaton and its decay products. In the considered model a stable compactification can be ensured via 
1 trapping of the modulus field by a minimum of the effective potential. 

JL PACS number(s): 04.50.+h, 98.80. Hw 

' 1 Introduction 



It is well known that the fundamental physical constants in superstring theories are related to the vacuum 
expectation values of the dilaton and moduli fields, and variations of these fields would result in variations of 
the fundamental constants. In the context of standard Kaluza-Klein models moduli are defined by the shape 
and size of the internal spaces (we shall refer to the corresponding fields as geometrical moduli) . Up to now, 
there are no experiments which show a variation of the fundamental constants (theoretical and experimental 
bounds on such variations can be found e.g. in Jl]]). According to observations the internal spaces should 
be static or nearly static at least from the time of recombination (in some papers arguments are given in 
favor of the assumption that a variation of the fundamental constants is absent from the time of primordial 
nucleosynthesis) . Therefore, part of any realistic multidimensional model should be a mechanism for moduli 
stabilization. 

Within superstring theories such a stabilization is achieved, e.g., via trapping of the moduli fields by 
Kahler [[| or racetrack 3| potentials. Recently it was pointed out in Ref. jl[], that in a cosmological setting 
the trapping stabilization of the dilaton field can enhance and become robust due to the coupling of the 
dilaton to the kinetic energy of ordinary matter fields. 

Within multidimensional cosmological models of Kaluza-Klein type the problem of geometrical moduli 
stabilization by effective potentials was subject of numerous investigations H. It was shown that the 
scale factors of the internal spaces can stabilize, e.g., in pure geometrical models with a bare cosmological 
constant and curved internal spaces as well as in models with ordinary matter. Small conformal excitations 
of the internal space metric near the minima of the effective potential have the form of massive scalar fields 
(gravitational excitons) [jq] developing in the external spacetime (later, since the sub-millimeter weak-scale 
compactification approach these geometrical moduli excitations are also known as radions). Gravexitons 
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were investigated for a number of models in Refs. J7|. In Refs. |j] - [J} the 4-dimensional Planck scale 
Mpi(4) was implicitly understood as the D— dimensional fundamental scaleQ and it was assumed that the 
internal spaces are compactified at sizes somewhere between the Planck scale Lpi ~ 10 _33 cm and the Fermi 
scale Lf ~ 10~ 17 cm to make them unobservable. 

Recently it has been realized that the higher-dimensional fundamental scale M Pl ^ D ^ can be lowered from 
the 4-dimensional Planck scale Mpi^ = 1.22 x 10 19 GeV down to the Standard Model electroweak scale 
Mpi($ + D') ~ Mew ~ 1 TeV providing by this way a new scenario for the resolution of the hierarchy problem. 
The corresponding proposal led to an intensive study of various multidimensional theories. The considered 
models can be roughly devided into two topological classes. 

The first class consists of models with warped products of Einstein spaces as internal spaces. For sim- 
plicity, the corresponding scale (warp) factors are usually assumed as depending only on the coordinates 
of the external spacetime. Whereas gravitational interactions in such models can freely propagate in all 
multidimensional (bulk) space, the Standard Model (SM) matter is localized on a 3T>rane with thickness of 
order of the Fermi length in the extra dimensions. Such a model was used in Ref. J8| for the demonstration 
of the basic features of the sub-millimeter weak-scale compactification hypothesis with Mpim) ~ Mew- 
Different aspects of geometrical moduli stabilization for such models were considered e.g. in Refs. |J. A 
comparison of effective cosmological constant and gravexciton masses arising in the electroweak fundamental 
scale approach with those in a corresponding Planck scale approach was given in Ref. juj]. 

The second class consists of models following from Hofava-Witten theory jllj where one starts from 
the strongly coupled regime of Eg x Eg heterotic string theory and interprets it as M-theory on an orbifold 
R 10 x S 1 1%% with a set of Eg gauge fields at each ten-dimensional orbifold fixed plane. After compactification 
on a Calabi-Yau three- fold and dimensional reduction one arrives at effective 5— dimensional solutions which 
describe a pair of parallel 3-branes with opposite tension, and location at the orbifold planes. For these 
models the 5— dimensional metric contains a 4— dimensional metric component multiplied by a warp factor 
which is a function of the additional dimension The geometrical moduli stabilization for such models 
was considered, e.g., in Refs. | JL3| . Difficulties of the stabilization in such models connected with the 
required unrealistic fine-tuning ofthe equation of state on our 3-brane were pointed out in Ref. 114], with 
possible resolutions proposed in ]l5fl . Moreover, in Refs. Jl6| ] it was shown that the stabilization of the extra 
dimension is a necessary condition for the correct transition from 5— dimensional models with branes to the 
standard 4— dimensional Friedmann cosmology. Various other aspects of cosmological brane world scenarios 
were investigated, e.g., in Refs. |17| . 

Usually, the moduli stabilization is based on the trapping of the geometrical moduli fields at a minimum 
of an effective potential so that the fields are static (or may at most oscillate near this minimum due to 
quantum fluctuations). However, in absence of any minima, nothing forbids them to evolve very slowly as 
long as their evolution does not contradict the observable data. It is clear that a sufficiently slow evolution 
is allowed, if the moduli fields came into this regime bef ore primordial nucleosynthesis. Such an approach 
is in spirit of Paul Steinhardt's quintessence scenario |l8|]. It may happen that moduli fields asymptotically 
tend to some limit. We shall call such a behavior dynamical stabilization. An example for a dynamical 
stabilization in modular cosmology was presented in Ref. |19| ]. However, the interaction with ordinary 
matter fields can destroy this stabilization mechanism. The situation occurs, e.g., in the model considered 
in sections ^| and |^ of the present paper. 

Subject of the investigation in the present paper is a multidimensional cosmological model with a minimal 
coupled scalar field as inflaton field. In Ref. |^l| it was pointed out that a dynamical stabilization of 
the geometrical modulus could be possible for a model with one Ricci-flat internal space and a zero bulk 
cosmological constant. Here, we investigate this model in more detail and show that the interaction of the 
geometrical modulus field with the inflaton field in most cases destroys the dynamical stabilization and leads 
to decompactification of the extra dimension. There are only two possible ways for a stable compactification 
of the extra dimension. For the first one, a dynamical stabilization, one has to assume that the modulus 
field is only coupled to the inflaton field but not to its decay products. Due to the exponential decay of 
the inflaton during reheating the force term in the modulus equation, obtained from the effective potential 
of the model, decreases also exponentially. The present friction term provides then the stabilization of the 
modulus field. The second possibility for a stabilization consists in the trapping of the modulus field near 
a minimum of the effective potential. As simple example we consider the effective potential from Ref. fllOl] 
where the minimum is generated by a non- Ricci-flat (curved) internal space and a non-zero bare cosmological 
constant. The analysis does not depend on the choice of the D— dimensional fundamental scale. 

The paper is organized as follows. In section ^ we explain the general setup of our model and show a 
possible mechanism for a dynamical stabilization of the geometrical modulus in zero-order approximation. 
In section we investigate the dynamical behavior of the modulus and the inflaton field in more detail and 
show that the interaction between them in most cases results in a decompactification of the internal space. 
An example for a stable compactification of the internal space by trapping of the modulus field near the 
minimum of the effective potential is presented in section kl The brief Conclusions of the paper (section 
B|) are followed by an Appendix on higher dimensional perfect fluid potentials and specific features of their 
dimensional reduction. 

1 Hereafter, D = D' +4 is the total dimension of the multidimensional spacetime, D 1 — the dimension of the internal components. 
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2 Model and general setup 



We consider a cosmological model with metric 

(0) . 23(x) (1) _ (0) . 7.2/ s (1) / n ,\ 

g = g y > + e py >g y > = g y > + b (x)g y > , (2.1) 
which is defined on a manifold with warped product topology]^] 

M = M x Mi , (2.2) 
where x are some coordinates of the Do = {do + 1) - dimensional manifold Mo and 

g m =g$(x)dx l, ®dx u . (2.3) 

(The model corresponds to the type one class discussed in the Introduction.) 

With D — Do + di as total dimension, n 2 D a D— dimensional gravitational constant and Sygh the 
standard York - Gibbons - Hawking boundary term we choose the action functional in the form 



S 



±- j d D x^/\g~\R[g] - 1 J d D x^/\g~\ (g MN Om^nX + 2£/(x)) + Sygh , (2.4) 



where the minimal coupled scalar field \ with an arbitrary potential U (x) depends on the external coordinates 
x only. This field can be understood as a zero mode of a bulk field. From the other hand, such a scalar field 
can naturally originate also in non-linear D— dimensional theories |^2j where the metric ansatz (2.1) ensures 
its dependence on the x coordinate only. 

The main goal of the present paper consists in an investigation of the internal space dynamical stabiliza- 
tion. As it was shown in |2l| ], such a possibility exists for the considered model only if a bare D— dimensional 

cosmological constant A identically equals zero and the internal space is a Ricci-flat one: R [ff^J — 0. 

Let 6o = Lpie 130 be the compactification scale of the internal space at the present time and 



v dl = Vq x Vi = bg 1 



d dl yj\g^\ (2.5) 



the corresponding total volume of the internal space. Instead of j3 it is convenient to introduce the shifted 
quantity: 

= 0-00. (2.6) 
Then, after dimensional reduction action (2A) reads 

S = i J d D °xJ\gV)\e d ^ {r[ 9 ^] -d l (l-d l )g^d»0d„0-g^ K ld„xd»X-2KlU(x)} 

M 

= / d D «x^\^R [ 5 C°>] - - Kl^g^d.xduX - 2*lmx)) , (2-7) 



where 



e dlP 

2kq 2k§ \bo J 

and we redefined the scalar field \ and its potential as follows 



b \ dl 



= (2-8) 



X = V^iX, V(x)=v dl U( X /V^) ■ (2-9) 



In action (2.7) n is the Do— dimensional (4-dimensional) gravitational constant: 



^ = T7T^— » ( 2 - 10 ) 

M Pl «d! 

where Mpi — Mp;( 4 j = 1.22 x 10 19 GeV. It is clear that the scale of the internal space compactification bo 
is defined now by the energetic scale of the D— dimensional gravitational constant k 2 d . If we normalize kJj 
in such a way that njj = Hir/M^y 1 , where Mew ~ ITeV is the SM electroweak scale, then we arrive at 



8tt 



This mprlplxan be easily generalized to manifolds with warped product topology of n internal spaces: M = Mo X Mi X ... X M n 
(see Refs. |§ |Uj). 



3 



the intensively discussed estimate bo ~ iQ( 32 / d i) "cm. (Here, we assume for a toroidal internal space 
Vi = 1.) In the present paper we shall not discuss different choices of k^. It has no influence on the main 
conclusions about the possibility of an internal space d ynar ni cal .s tabilization . Generally speaking, we might 
not fix a scale of the internal space stabilization in Eqs. (|2.6[ ) - fl2.8|) , but simply set there /3o — and 60 = Lpi . 
For our concret model this does not mean that the stabilization takes place at the planckian scale but results 
simply in a planckian normalization of the D— dimensional gravitational constant: n 2 D ~ Mp^ 2+dl ^ jnl. In 
the latter approach, the scale of stabilization is defined from the equations of motion. 

To get the equations of motion in the Brans-Dicke frame it is convenient to rewrite action (p. 7]) as follows: 



M 



(2.11) 



where cj := (1 — di)/di < and 



(0) 



Varying this action with respect to the metric g^J and the fields $ and x, we get respectively 



and 
where 



z? (0) e> _ V _i_ w 



n x +±g^^d fl ^x = V , 



(2.12) 

(2.13) 
(2.14) 
(2.15) 



T 

a = 



T£ = D L m +g m ^d^ X d„X, 
=9^^\g^^ 



V\g 

and the prime d enote s the derivative with respect to the field x- 
Contracting (2.1J) with respect to g^ M " yields 



(2-D )/2 J + <VW + $(D - 2 )/2 UV - 



(2.16) 
(2.17) 

(2.18) 



Now, combining (2.14) and (2.18) we can rewrite the equation of motion for the field $ in the convenient 
form 



2oj + 



Do - 1 
(D - 2)/2 



T — 2KoL m 



D -2 



Do -2 



(2.19) 



Let us specify metric g^ as 

5 (0) = -dt®dt + a 2 (t)g {0) , (2.20) 

where g*- ' is the metric of the do— dimensional Ricci-flat external space: R ^ < "°' ) j — kdo(do — 1) = => fc = 
(in accordance with recent observations our Universe is flat). For this metric we get 



Roo[g 



(0)1 



-do- and R[g {0) ] = do 



2- + (do -l)-5 



(2.21) 



where the dot denotes the derivative with respect to the synchronous time t in the Brans-Dicke (BD) frame 
(string frame). 

Below we consider homogeneous fields: $ = $(i) and x = xM- Then, the 00— component of the Einstein 
eq. ( 2.1S ) reads 



H 



2 $ 
-H 



d -l $ d (d -l)V*y do(d -l)V2 



+ 



2^o 



(2.22) 
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where H = a/ a is the Hubble parameter. If we take into account th at for our metr ic (2.2C) and any 
homogeneous field <f> holds 0<f> = — <f>— doHcj), then we get for equations (2.19) and (2.15) correspondingly: 



$ + d H$ = 



2k, 2 , 



uj(d - 1) + d 



and 



X - 



d Hx +$X = -V . 



(2.23) 



(2.24) 



We suppose that the potential V(x) of our model has a zero minimum: V\ min = 0. Then, after the \ 
field has evoluted down-hill to this minimum and is frozen out, the right hand side of eq. ( |2.23| ) is equal to 
zero and the equation has a solution 



*(«) = 23 - *> 



— d J Hdt 

dte 4 ° 



(2.25) 



where $o := $(t = to) is an initial value for $. The constant of integration in ( 2.25 ) is chosen in such 
a way that <3>(t — ► oo) 



l/2«o, i-e. that it corresponds to the stabilization of internal space b 



Thus, in zero-order approximation (with respect to the interaction between the $ an d y fi elds in eq. (2.23)) 
the dynamical stabilization of the internal space takes place if the integral in eq. (2.25) is convergent at 
the upper limit t — > oo. If the integral diverges the dynamical stabilization mechanism fails to work and 
decompactification occurs. A detailed analysis of this problem will be given in the next section where we use a 
self-consistent approach and investigate the influence of the interaction between the fields on the dynamical 
stabilization of the modulus field It is more convenient to perform the corresponding analysis in the 
Einstein frame. Obviously, if stabilization occurs in the Einstein frame it has place also in the Brans-Dicke 
frame and vice versa. 



3 The Einstein frame 



A conformal transformation to the Einstein frame is given by equation 



(3.1) 



and yields 



where 



M 



<P ■ 



«0 



1 Jdi{D-2)- 



D -2 



and the effective potential can be written as 

U e ff[<p,x] ■= e 2K ° V V 



(3.3) 



(3.4) 



(For definiteness we use the minus sig n in eq. (3.3).) 

The dimensionally reduced action ( |3.2| ) describes a system of 4-dimensional gravitational and scalar fields. 
The problem of the internal space stabilization is reduced now to the investigation of the dynamics of these 
fields. For this purpose we specify the metric g' '. In the Einstein frame the 4-dimensional metric ( |2.2C| ) 
can be rewritten as 

g {0) = -dt®di+a 2 (t)g {0) , (3.5) 
where 'tilded' quantities are related to the Einstein frame^. 



Then, the equations of motion following from action (3^) read: 



H 2 



2k 



t;^ 2 + + u *ff 



(p + d Hip 



do(do-l) I2 r '2 



dtp 



(3.6) 
(3.7) 



5 



X + d Hx 



ft 



dx 



(3.8) 



Hereafter, dots denote derivatives with respect to the synchronous time t in the Einstein frame, and primes 
denote derivatives with respect to the infiaton field x- Further on, we assume for the external space dimension 
the usual value do = 3. For this choice of do the parameter a is fixed by the dimension of the extra space as 
Vrfi/2(di + 2). 



Inflation in our particular model (3.6) - ( p.q ) is known as soft inflation and was extensively studied in 
Refs. |2(], Ell. Here we simply quote the corresponding results. During inflation the slow roll conditions 
ip 2 < 2U ef f, x 2 < 2[/e//, (p < 3Hip and x < 3Hx hold and the Eqs. (3~6) - ( p3~g| ) reduce to the simpler 
form 



SHip 
3H X 

H 2 



-2<TKoe 



V{X): 



2 

K-0 



u, 



v'(x), 

2 

^£ e 2CTK0V l / (x) ■ 



(3.9) 



As result, the fields x an d V? are connected with the scale factor of the external space by the relations 



—H 

Kq 



Ko(<p — ipo) ~ — 2ct / Hdt = —2aln — 

do 



(3.10) 



and 



2 
«-() 



xo 



0(1 



(3.11) 



where the initial conditions a re d en oted by the subscript " 0" and are fixed at to , when the Universe enters 
the inflationary phase. From (p. 3), ( 3.1C ) we see that during the inflationary stage of the external spacetime 
the scale factor of the internal space undergoes inflation too (see e.g. also 



bit) 



u t^-tt ! Hdt 

One - 2 J 



(3.12) 



The necessity of such a stage was stressed in paper |23| , where it was shown that to solve the horizon and 
flatness problem, there must be a stage of inflation in the bulk space before the compactification of the 
internal spaces can be completed. If we take into account that during this stage a ~ exp ( f Hdt) then we 
get following relation between initial and final values: bf/bi = [aj /a,i\ 2 ^ D ~ 2 ^ (see also pl|). 

Let us investigate now the dynamical behavior of the model at the stage when the infiaton scalar field 
X is evoluting down-hill to the minimum of its effective potential and is located in the very vicinity of this 
minimum. Without loss of generality, we suppose that potential V(x) has a zero minimum at the point 
X = : V\ x=0 = 0, V'\ x=0 = , V"\ x=0 > 0. This leads to two implications. 

First, for |x| < 1 (i.e. jxi < Mpi) this potential can be crudely approximated by 



where m. 



V"L_ . In the same approximation we get for eqs. (3.7) and (3. £) correspondingly: 
ip + 3Hip + KQae 2<TKOV 'm 2 : x 2 = 



(3.13) 



(3.14) 



and 



X + 3#x + e 



av m x x- 



0. 



Second, from the structure of the effective potential UeffiXy 1 ^) 



(3.15) 

e 2<TK ° v 1/(x) it is clear that the zero 



minimum V\ _ = is globally degenerate and, in crude analogy with Goldstone bosons in ip —theories, 



the field tp plays the role of the zero mode along the degeneration line x 
easy to see from the Hessian of the effective potential 



in the (x, <p)— plane. This is 



d 2 

4 



Ueff\ 



X = 





e 2aK0V m 2 x 



(3.16) 



which defines the mass matrix of the normal modes of the model at the point x — (see also Ref. |6J). This 
means that the geometrical modulus field tp corresponds to a flat direction of the effective potential U e ff and 
is not stabilized by a minimum of this potential (see e.g. Ref. j24|). Due to the remaining ip— dependence 
of the effective mass of the x~ mode the analogy with Goldstone bosons in ip 4 — theories is only very crude. 



() 



As first step of our analysis, we in vest igate the equation of motion (3.14) of the modulus field not taking 
into account the energy constraint (3.6) and assuming that the \ field is frozen in the minimum of the 



potential V(x) (zero-order approximation in the fluctuations of field \) 



X 







(p + 3Hip « . 



(3.17) 



For a cosmologically non-damped evolution with H = we would have ip as and a dynamically non- 
stabilized behavior of the modulus field ip = ipot + ipo, where the initial val ues a re chosen at some i = to : 
ipo := (p(i = to), ipo := <p(t = to). Due to the damping term 3Hip ^ eq. (3.17) has a solution 



T -3 / Hdt 

■pit) = -ipo I die { o 
t 



(3.18) 



which describes a dynamical (as ympt otical) stabilization if the integral converges at its upper limit t — ► oo. 
The constant of integration in (3.15) is chosen in such a way that ip(t — > oo) — > and corresponds to 
the internal space stabilizati on at b — > bo- Obviously, this solution is the Einstein frame analogue of the 
Brans-Dicke frame solution ( |2.25| ) and reproduces the dynamical stabilization scenario of Ref. |HJ. For a 
power-law behavior of the external scale factor: a ~ t s => H = a/a = s/t (s = 1/2 , 2/3 corresponds to the 
radiation dominated (RD) and matter dominated (MD) era respectively) this equation yields 



<p(t) 



1 3s 



3s - 1 t : 



1 

S> 3 



(3.19) 



We can use this equation for estimates of the variation of the fundamental constant s due to the dynamics 
of the internal space. Usually jjj, such variations are proportional to b/b. From eq. ( 3.19) we get 



D 



d 1 (D-2f° ip \t 



to 



(3.20) 



We choos e the initial values ipo and to in suc h way that they correspond to the end of inflation: Koipo = —2aH e 
(see eq. and H e t ~ 70. Then, eq. (|3.2cj) yields 



D- 2 



70 

Het 



(3.21) 



For the simplest models of inflation, as in our case, the COBE data predict H e ~ 10 °Mpi. Thus, we 



obtain b/b\ _ 



10" 14 yr _1 - This means that effectively there is no variation of the fundamental 



li~10 3 sec;s = l/2 

constants starting from the time of nucleosynthesis. For the modulus field value we can easily get an estimate: 



Kotp(i) 



140(7 / - A 3s — 1 
'70/ffet) 



3.s 



3 = 1/2 



10" 18 < 1. 



(3.22) 



Thus, in zero order approximation and not taking account of the constraint (3.6), the internal space would 
appearently stabilize to this time. 



Let us now use a self-consistent approach to our model taking also into account the constraint (3.6). For 
frozen \ ne ld we have x — Xo = 0, \o = and V( yo) = so that the energy density P x (e) = \x + U e ff 
vanishes {p x (E){Xo) = 0) and the equation system (3.6) - (3.8) reads: 



The solutions are easily found as 
1 



H = 



rl/3 



3t 



ip + 3Hip 
nop — ± 



2 

K0 -2 

= 0. 



2 1 

3? : 



Kop 



2 , t 

- In j- + K ipo 

3 to 



(3.23) 
(3.24) 

(3.25) 



and show that the possible dynamical stabilization mechanism, which holds for s > 1/3, fails to work 
in the case s = 1/3 due to the logarithmically diverging ip. Heuristically, we can say that for frozen 
inflaton field x the decoupled modulus field ip behaves like an ultra-stiff perfect fluicj^ with equation of state 

See ]25[ |26| and the Appendix for a discussion of the equivalence between a scalar field ip and a perfect fluid with energy density 
p v , pressure iP v and equation of state P w = (a — 
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Ptp(E) = Pif>(E) ~ 2 6 an d produces not enough " cosmological friction" H — s/t'm order to come to a rest at 
some finite value \<p(t — > oo)| < oo. As consequence the additional dimensions cannot stabilize, but rather 
they decompactify. 

This behavior can be partially circumvented by passing from the homogeneous modulus approach <p = 
p{t) to an inhomogeneous one <p = <p(x) with modulus fluctuations behaving like radiation (see e.g. jl9|, Bj|). 
The corresponding energy density could lead to the needed "cosmological friction" H with s > f/3. But as 
we will show below in this section, the interaction of the modulus field with the inflaton will destroy such a 
dynamical stabilization mechanism. For simplicity, we will restrict our subsequent considerations to purely 
homogeneous fields 



For such fields we can circumvent the decompactification mechanism (3.25) if we assume that x performs 
small fluctuations around the minimum of its potential V(\) yielding by this way a nonvanishing energy 
density P x (e) > which could provide the nee ded " cosmo lo gical friction" H with s > 1/3. So, we shall 
investigate corrections to the equation system fl3.6| ), (3.14), ( |3.1q ) due to the dynamics of the field \. To 
achieve this g o al, we assume that the modulus field ip is already nearly stabilized at ip « and embed our 
system (3.14), (3.15) in a generalized astrophysical setting allowing for decay processes of the inflaton field 
X into usual matter. The corresponding part icle interactions result in a polarization operator II of the field 
X which shifts the squared mass in eq. ( |3.15| ): m x — > m x + II (see e.g. |28|). The imaginary part of II is 
responsible for the inflaton decay and is connected with the decay rate T x oy the relation Imll = m x F x . It 
is shown in Ref. |^8| that phenomenologically such a dec ay can be taken into account by adding an extra 
friction term Tx to the classical equation of motion ( |3.15| ) (instead of adding the term proportional to the 
imaginary part of the polarization operator): 



X+(3H + F x )x + e 2 ^m 2 xX = 0. 



(3.26) 



We use this equation to crudely understand the dynamical behavior of the field x during the post-inflationary 
evolution of the Universe (not taki ng into account preheating or subtleties of the decay processes pq| ). It is 
convenient to investigate eq. ( |3.26| ) with the help of a substitution (see J2!| , Chapter 14) : 



X (t) := B(t)u(t) := e 



FOV+3/ Hdt) 



U(t) . 



(3.27) 



where the function u(t) satisfies the equation 



u + 



rh 2 x --(3H + F x ) 2 2 



u = 



(3.28) 



and rntit) = m x exp(2aK,op). If we suppose further that (3H + F x ) 2 , H, (m x /m x ) 2 are small compared 
with m?. (which, for a viable model should take place in the large t 3> io limit when \nop\ <C 1), then eq. 



(3.29) 



( p.28|) has an approximate solution of the form (see also [ [jQfl ) 

u(t) ~ cos(rh x t) x(t) ~ B(t) cos(fh x t) . 

It can be easily seen from the definition of B(t) that this function satisfies the equation 

d /~3r}2\ -3 D 2 

— [a B = — 1 x a B . 
Approximating the energy density of the inflaton and the corresponding number density as 



(3.30) 



Px(E) = \x + \m 2 xX 2 
shows that they satisfy the differential relations 



1 D 2 -2 

— B 



L X(E) 



— B m-v 



(3.31) 



?( a Px(E)) = ~r x a p x(E) 



and 



-Ja n x(E) ) 



= —1 v a n 



X(E) 



(3.32) 



with solutions P x (b) ~ e~ Txt a~ 3 and n x (E) ~ e~ Fxt a~ 3 . From these relations we see that during the stage 
fh x > H 3> F x the inflaton performs damped oscillations with an energy density corresponding to a dust-like 

perfect fluid {p x (E) ~ 2 -3 ) with slow decay ~ e~ Txt ~ 1. This reheating stage ends when H < F x and the 
evolution of the energy density is dominated by the exponential decrease due to decay with rate F x . 

If the decay rate F x of the inflaton particles into usual Standard Matter (SM) is sufficiently large, then 
starting from the characteristic decay time to ~ T" 1 , i.e. from the time of the most intensive reheating, 
the energy density of the corresponding relativistic particles behaves as psM(E) ~ 5-~ 4 . Clearly, the energy 
loss of the x field due to the decay process is accompanied by a corresponding energy increase of the decay 
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products. As result, the effective energy density p e //(E) := P x (e)+Psm{e) is only deluted by the cosmological 
expansion and can be roughly approximated as 



Peff(E) - P X (E) + PSM(E) 



<*a- 3 + (e" 



(3.33) 



This relation holds for times t ^ ti ~ ^ x with <§C 1, where ti plays the role of an effective initial 

time which fixes the begin of the coherent \ oscillations and of the decay process. From the Einstein 
equations of the extended interacting system we can derive a corresponding Friedman equation similar to 
(3.6). Obviously, for a viable model the internal space should be stabilized before the nucleosynthesis, i.e. 



parameter is defined from this extended Friedman equation by p e jf( E y. H 2 ~ ^KoPef f(E)- Thus, the Hubble 
parameter during the post-inflationary stage, including the period of nucleosynthesis, is defined by the 
energy densities p x (E) and Psm(e) and can be roughly approximated as H = s/t where s = 1/2,2/3 for 
p e ff(E) ~ 5 -4 , 2 -3 respectivelyQ Hence, in the same approximation we have 



Peff(E) « —j-t 



(3.34) 



Let us now in some detail consider the required asymptotic evolution of the modulus field ip ~ during 
the post-inflationar y stag e up to the period of nucleosynthesis, when the modulus stabilization should be 
finished. From eq. ( A. 13 ) give n in t he Appendix we see that during the inflaton dominated post-inflationary 
stage the equation of motion j3.14 ) for the modulus field ip can be rewritten as: 



Cp + 3Htp = (a x - 2)an p x (E) 



(3.35) 



After this stage, the further evolutional behavior of the modulus field crucially depends on its coupling to 
the decay products of the inflaton field, i.e. on the coupling type to the SM fields. We will illustrate this 
specific model dependent feature with the help of the two most simplest examples, a model (I) where the 
decay products of the inflaton have a similar coupling to the modulus field like the inflaton itself, and a 
model (II) where the decay products are not coupled to the modulus field at all. 

Model (I): Under a model with similar coupling of the inflaton field and its decay products (the SM 
fields) to the modulus field we understand a model where the evolution of the modulus field is defined via 
extension of (3.7) by an equation of motion of the type 



Cp + 3Hp = (a x — 2)aKop x (E) + (asm — 2)okopsm(e) 

~ (a e ff-2)<TKop B ff(]ny 



(3.36) 



Let us furth er ass ume that the required asymptotic stabilization of the modulus field is almost achieved ip w 
so that eq. (3.36) c an he considered as a non-homogeneous differential equation. In this approxi mation the 
3) can be written as a sum of the general solution of the homogeneous eq. (|3.17[) 



and 



general solution of ( p. 36 

a particul ar so lution of the n on-ho moge neous equation. As a solut ion q f the homogeneous equation we can 
take eq. ( |3.19| ). Then, using (3.34) and (A. 17) the solution of eq. ( |3.36 ) can be approximated as 



<po to* 
3s-li3 3 - 



2so 

HQ 



111 



t 



(3.37) 



with s = 1/2, 2/3 for p e ff ~ a 4 ,a 3 respectively. The time To plays to role of an effective initial time which 
is d efined from the val ue of the internal scale factor b(t) at the beginning of the evolutional stage described 
by ( |3.36| ). From eq. ( ]3.37j ) we see that the used perfect fluid ansatz leads to a decompactification of the 
internal space. For times t ^> to, the internal scale factor behaves as b w bo (t/To) 1 , where 7 = 2s/(D - 2) 
and e.g. 7 = 1/4,1/3 for d\ = 2 and s — 1/2,2/3 respectively. Eq. (3.37) shows also that for t ^> to 



we obtain b/b ~ —Ka<p ~ i , so that via extrapolation to our present time we would get the esti mate 
10 _10 yr _1 . This estimate is much greater than 10 _14 yr~ 1 following from observations. So, solution (3.37) 
leads to a decompactification of the internal space and does not ensure the necessary stabilization of the 
internal spaces at the present time. This was proved by the rule of co ntraries. First, we supposed that the 
modulus stabilization occurs and we can use the approximation (3.34). Then, we showed that our proposal 
is wrong because an unboundedly increasing destabilization term occurs in the solution (3.37). 

Let us now consider model (II) w ith vanishing coupling of the decay products of the inflaton to the 
modulus field^]. In this case eq. (3.35) of the modulus evolution holds also after the inflaton-dominated era. 



4 We note that the expansion parameter s is connected with the parameter a in the equation of state P e ff(E) 
by the relation s = tt—- See also (A. 15) - (A. 17). 



(a - l)p e //(B) 



3 But the decay products still define the dynamical behavior of the Universe. 
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Using the approximation (3.34) for the energy density p x o f the decaying dust-like perfect fluid component 
we rewrite this equation with an ansatz analogous to (|3 .34) and s — 2/3 as 



ifi + 2i~ 1 ip= -p-t~ 2 e - T ^. (3.38) 



For simplicity of notations we introduce the abbreviation q := 4cr/(3Kor x ). Then, the solution of (3.38) with 
internal space stabilization tp(t — > oo) — > at & — > bo and initial condition ip(t — to) = tpo can be easily 
found as 



tpit) = [qe- T ^ - tpoPo) J- 9 J. dt • ( 3 - 39 ) 

The corresponding characte ristic variation of the internal scale factor for the same initial conditions on to 



and ipo as for the estimate (3.21) reads 



and gives for a decay channel of the inflaton to fermions with decay rate F x ~ 10~ 12 Mpi the esti- 
mate &/& ~ 10 _31 j/r _1 . The modulus field stabilizes at the same time at Ka<p\j ln 3 r , ~ 

U~~10 3 sec,di=2 sec, ai 

10 -35 <C 1. Thus, there exists a possible dynamical stabilization scenario for a decaying inflaton field and 
a modulus field which is not coupled to the decay products. Clearly, this dynamical stabilization scenario 
is rather artificial and in general the decay products will be also functionally coupled to the modulus field 
what can destroy the stabilization. 

Above, we considered a model with zero effective cosmological constant. However, resent observations 
show the existence of a positive cosmological constant A ~ lCP 57 cm~ 2 for our Universe. So, it is of interest 
to include such a A— term into our consideration. This can be easily done if we suppose that the potential 
V(x) has a non-zero minimum at x = 0. Thus, for \x\ J$ Mpi the potential reads 



V(X) ~Vo + \m 2 xX 2 • (3.41) 



It is clear that A e ff = e 2<TKO<p Vo plays the role of the effective cosmological constant which asymptotically 
tends to Vo in the case of the internal space stabilization if — > 0. Such a behavior of the effective cosmological 
constant is similar to the quintessence scenario flq]. We suppose that in accordance with observations 
KqVo ~ 10 -57 cm -2 > 0. For such values of Vo the influence of Vo on the Universe and the fields dynamics 
becomes essential only at a stage close to our present time. During earlier post-inflationary evolution stages 
e 2tTK ° v Vo is negligible compared with P x (e)'- e 2tTK °' f Vo P x (e) = |x 2 + ^ r "-xX 2 - But, with progression of the 
Universe expansion P x (e) decreases and becomes less and less compared with e 2tTKOV Vo. Let us investigate 
the influence of Vo on the internal space stabiliza tion when e 2tTK ° v Vo > P x (e)- We also assume that up to 



this time |ko<p| <C 1. Under this assumptions eq. (3.35) is modified as follows: 

(p + 3H(p = -2Koae 2cTKW Vo + (a - 2)n ap x{E) » -2k ctVo . (3.42) 



A solution of this equation can be found similar to eq. (3.37) and reads 



= ^_ JiL _ . (3.43) 

^ v ' 3s -It 3 '- 1 3s + 1 K ' 

Thus, an effective cosmological constant results also in a decompactification (destabilization) of the internal 
space at late times. 

Summarizing the results of the present section, we can conclude that the considered possible dynamical 
stabilization scenario is very sensitive to the coupling of the modulus field to small fluctuations of the inflaton 
and matter fields, as well as to a non-vanishing vacuum contribution. Thus, this scenario is in general not 
sufficiently robust and each concret model needs a detailed study whether a dynamical modulus stabilization 
could work or not. In the next section we extend our toy model and induce a trapping mechanism to guarantee 
the modulus stabilization. 



4 Stable compact ificat ion 



In this section we present an example for a model which ensures stable compactification of the internal 
space by a trapping of the geometrical modulus at the minimum of the effective potential. For this purpose, 
we modify our setup model of section H including a non-zero bare D— dimensional cosmological constant 
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A into the corresponding action functional 
R[g w ] = -Ri = const ^ 0. Then, instead of 



(2.4), (2.7) and assume that the internal space is not Ricci-flat 
(3.4) we get an effective potential 



Ueff[<P,x] 



= 2 "W\/ (D-2HD -2) 



d 1 (D-2) 



+ -,A + V(x) 



2<TK tp 



2kq 



Rie 



2~fK ip 



7 A + V( X ) 



(4.1) 



where Ri := Ri/bg and we used the obvious abbreviations a = \J d\j2(d\ + 2) and 7 = y/2/di(di + 2). 
We restrict ourselves to the case of a potential V(x) with a single zero minimum V(xmin) = V(x)\ ■ = 0, 
dxV\ min = m x > because the non-zero minimum case V(x)\ min = Vo is trivially reduced to the zero 
minimum one by inclusion of Vo into A/kq- 

It is easy to check that the effective potential ( |l.l| ) has a global minimum at the point ip = 0, y = Xmin 
if the bare cosmological constant and the scalar curvature of the internal space are both negativerj A < 0, 
Ri < 0. Additionally it is necessary that these parameters are connected with the compactification scale foo 
by a fine-tuning condition 

„ . . . -Ri , D — 2 R\ 

2aA = (a + 7 )- r • ' 



b 2 o 



A = 



2di 



I, 2 
°o 



From the Hessian (the mass matrix) of the effective potential 

St 



d 2 



U, 



eff\ v= o 7Xv 



2 7 (o- + 7)i?i 




(4.2) 



(4.3) 



it is clear that in contrast with ( p.lfiD the minimum is non-degenerate and induces a non-vanishing mass 
m 2 , = — 2^(a + r y)Ri of the modulus fieldf^]. As result the position ip = is energetically favored and provides 
the necessary trapping. 

The trappin g ca n also be seen directly fr om the equations of motion of the (f,x) syste m. W hereas 



substitution of (4.1) into the x field equation (3 
modulus field reads now 



does not change this equation, the equation (3.7) for the 



ip + doHip 



dU, 



eff 



dip 



(cr + 7).Rie 



^A- 

K0 



2aV( X ) 



(4.4) 



So, for fine tuned parameters (4.2) the point 93 = 0, % = x-min is the trivial solution of the system (3.8), 
(kL4) and small linear perturbations around this solution satisfying 



Sip + doHSip — 2j(a + ^RiSip 



(4.5) 



6x + d HSx + m x &x = 

are damped by the present friction terms. The global modulus dynamics is easily understood from the 
form of the effective potential. The minimum at (tp = 0, x = Xmin) is the deepest point and defines after 
stabilization of the system the effective cosmological constant A e // = K^U e f /| . = Ri/di < of the exter- 
nal spacetime. In the small-scale-factor region we have U e ff(<p — > 00) — > 00 and in the decompactification 
region U e ffUp — > —00) — > 0, so that we are led to the following implications. 1) the Brustein-Steinhardt 
|33[ problemn cannot occur because there is no barrier in the model at all, separating the minimum of the 
potential from an asymptotically allowed region (the single minimum is located en erge tically deeper than 
the decompactification region); 2) higher order approximations to the solution of (4.4) will not alter the 
stabilization picture for the same reason; and 3) the dying away friction term doH<p> will leave at most small 
fluctuations of the geometrical modulus field around the minimum. This means, that this minimum stably 
compactifies (traps) the internal space. Clearly, such a stabilization mechanism is different from the above 
discussed dynamical stabilization where the effective potential has no minimum with respect to (p. 

As conclusion of th is se ction we would like to note that the negative effective cosmological constant 
A e // < of the model |y]) leads to a turning point in the scale factor evolution of the external space, i.e. 
a stage of external space inflation is followed be a period of contraction. This happens when the effective 
energy density p e // of the fields in the Universe becomes smaller than |A e //|. Model W^Ti) with stable 
internal space can be brought in agreement with the observed positive effective cosmological constant of the 
Universe, e.g. if one takes into account additional higher dimensional form fields p4|. 



6 Compact internal spaces with negative curvature R\ = —di(di — 1) can be constructed, e.g., as hyperbolic coset manifolds 
H dx/Y<l\ (for details see §2)). Here H d i is an infinite hyperbolic space and T dl — an appropriate group of discrete isometries. 
The coset manifold itself can be imagined to be built up from a fundamental polyhedron in H dl with faces pairwise identified. 

7 These massive modes (gravitational excitons HQ propagate in the external spacetime and can yield a considerable contribution 
to the dark matter in our Universe llOl . 

8 The Brustein-Steinhardt problem can occur for potentials with an energetically allowed decompactification region which is 
separated from a narrow valley around the minimum by a barrier. Supposing that the modulus field rolls down the potential it can 
overshoot this narrow valley and enter the decompactification region. The problem consists in the fact that the modulus field with 
high kinetic energy will not "feel" the valley and will not stabilize in the corresponding minimum. 
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5 Conclusion 



In the present paper we investigated in detail the recently proposed pifl mechanism of internal space dynam- 
ical compactification in a multidimensional model with one Ricci-fiat internal space and a massive scalar 
field acting as inflaton. As shown in section ^, such a dynamical stabilization could be possible if the mod- 
ulus field is only coupled to the inflaton but not to its decay products. Due to the exponential decay of 
the inflaton during reheating the force term in the modulus equation, obtained from the effective potential, 
decreases also exponentially so that the modulus field can dynamically stabilize via friction. If the decay 
products of the inflaton will couple similar to the modulus field like the inflaton itself, the dynamical stabi- 
lization will not occur. An analogous decompactification has place for a model with non-vanishing effective 
cosmological constant. A stable compactification for the considered model can be ensured via trapping of 
the internal space scale factor by a minimum of the effective potential. An example of such type of stable 
compactification is given in section ^. Oscillations of the modulus and the inflaton field around the minimum 
are observed as massive scalar particles (gravexitons and inflaton particles) in the external spacetime. Such 
particles were considered in the electroweak and the planckian fundamental scale approaches in Ref. [p"o||. 
In the same paper it was pointed out that there exists a rescaling for gravexciton masses in the different 
approaches and extremely light particles may arise in the planckian fundamental scale approach. Further it 
was shown there that particles with masses m v ~ 10 _33 eV are of special interest because, first, they do not 
overdose the Universe, second, the period of their oscillations T oac ~ l/m v is of order of the Universe age 
~ 10 18 sec, and third, the effective cosmological term A v = K^m^ip 2 ~ 10 -57 cm -2 corresponding to models 
with such particles and tp ~ Mpi, takes a value of order of the presently observable cosmological constant 
in the Universe. Thus, these particles evolve extremely slowly within their minimum and it is attractive to 
treat A v as an effective cosmological constant in the spirit of quintessence (for models with zero minimum 
of an effective potential). However, such a quintessence would have a drawback because it requires a highly 
fine tuned initial condition of the modulus field ip (see e.g. p5[|). 
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A The effective potential of a higher dimensional perfect 
fluid 

In this Appendix we assume that the bulk scalar field \ behaves like a perfect fluid living in the D— dimensional 
bulk spacetime and obeying the equation of state 



Px_ = (a- l)pz 



(A.l) 



with < a < 2 a constant. The energy de nsity p x = — T® of the scalar field x an d its pressure P x = 



are defined from the action functional (2.4). Further, we assume that the metric of the Dq— dimensional 



external spacetime is given by the ansatz{Z.20) and that a homogeneous a ppr oximation for the scalar fields 
X = xW) ^ — can De used. Then, as it was shown in |25[ , the action ( p.4| ) of the gravity — scalar field 
system is equivalent to the gravity — perfect fluid action 



S 



1 

2 K i) jm 



d D x x /\g~\\R[g]-2^p x (a,b)\ , 



where 



p x (a,b) = 



.4 



A := const 



and we took into account t he r edefinitions (2.9) of the scalar field \ — > x an d the potenti al U( y) 
Dimensional reduction of (A.2) gives the Brans-Dicke frame action functional similar to (2.11) 



(A.2) 
(A.3) 

v(x)- 



S : 



= J d D °x^\^~\{<i>R[g^] -^g^d^d^^24<fp x( B D) (a,<f)} , (A.4) 

M 



12 



where the energy density p x ibd) is defined as 



Finally, via conformal transformation (3.1) we pass to the Einstein frame and arrive at 



\ ( d Do xJwV\\\h\g™\ -~g^»d li <pd v <p-2p x{B) {a, V ) 



The energy density p x ( B )(2, 95) can be recast as 

Px(E)(a,ip) = 



o 2ctk 0V 



(2-a)o-K 45 



Px(SD)( a ' $ ) 

Px(a) , 



(A.5) 
(A.6) 

(A.7) 

(A.8) 

(A.9) 
(A.10) 



where we used the relations between the synchronous times dt — ±e aKoip dt and the scale factors a = e ""^o^q 
in the Brans-Dicke frame (t, a) and the Einstein frame (t, a) and defined the reduced energy density p x as 

.4 



Px(a) 



l ad obn 



(A.ll) 



The equations of motion following from action functional (A.7) for homogeneous field (p = ip(t) and 
energy density p x (E) = Px(E)(t)j an d Ricci-flat external space iJfg' -'] = read: 



tf 2 



2 K g 



+ Px(£) («>¥>) 



</5 + doHtp 
H 



do(do-l) 1 2 

= (a-2)ffK0P x ( B )(o,93) ; 



'Ml 



do 



- +ap x(J3) (a, <^)} 



(A.12) 
(A.13) 
(A.14) 



Eqs. (A. 10) and (A.13) show that for a = 2 the interaction between the modulus field tp and the perfect 
fluid (scalar field xjls absent. This fact is easily explained because a = 2 corresponds to a perfect fluid with 
ultra-stiff equation of state, which describes a scalar field with vanishing potential energy V(\) = 0. The 
other extremal case with a — describes a vacuum equation of state and can be used for potentials V(x) 
with non-zero minimum Vo- 

At the end of this Appendix let us consider a model with stabilizing modulus field tp ss and decaying field 

space. 




H 
H 



2^1 



do(d -l) PSM(B)(a) ' 



do - 1 

so that the Hubble parameter behaves as H = s/t with 



Psm{e)(o) ■ 



do a 



(A.15) 
(A.16) 

(A.17) 
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